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Abstract 

Synchronizability of stable, output-coupled, identical, time-varying 
linear systems is studied. It is shown that if the observability grammian 
satisfies a persistence of excitation condition, then there exists a bounded, 
time- varying linear feedback law that yields exponential synchronization 
for all fixed, asymmetrical interconnections with connected graphs. Also, 
a weaker condition on the grammian is given for asymptotic synchro- 
nization. No assumption is made on the strength of coupling. Moreover, 
related to the main problem, a particular array of output-coupled systems 
that is pertinent to much-studied consensus problems is investigated. In 
this array, the individual systems are integrators with identical, time- 
varying, symmetric positive semi-definite output matrices. Trajectories 
of this array are shown to stay bounded using a time-invariant, quadratic 
Lyapunov function. Also, sufficient conditions on output matrix for syn- 
chronization are provided. All of the results in the paper are generated 
for both continuous time and discrete time. 



1 Introduction 

When do the trajectories of a number of coupled individual systems converge to 
each other? This question outlines the multifaceted problem of synchronization 
stability. Unknotting this problem requires understanding the interplay of two 
pieces: the set of individual systems' dynamics and the (varying) topology of 
their coupling, i.e. who influences whom and how strongly. The general prob- 
lem is insuperably difficult, which has led people to a number of simplifications, 
justifiable for certain applications. For instance, when the individual system dy- 
namics is taken to be an integrator, by using convexity arguments, trajectories 
have been shown to converge to a fixed point in space as long as the (directed, 
time-varying) interconnection satisfies a fairly weak connectedness condition. 
Since, once synchronized, the righthand sides of the systems vanish, the word 
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consensus is used when referring to this case; see, for instance, [§1 EJ E H] • An- 
other direction of investigation is fueled by the fact that the speed/occurrence 
of synchronization is related to the coupling (strength) between the individual 
systems. Studies concentrated on understanding this relation have been fruitful 
and significant results have emerged. We now know that the spectrum of the 
interconnection matrix is where we have to look at if we want to measure the 
strength of coupling in order to determine whether synchronization will take 
place or not. Roughly speaking, under the assumption that some Lyapunov 
function (related to the individual system dynamics only) exists, one can guar- 
antee stability of synchronization if the coupling strength is larger than some 
threshold; see, for instance, [TH1 EH HE1 12 • There are numerous other interesting 
research directions accommodating notable works in synchronization stability. 
We refer the interested reader to the surveys [TH [17], [4[ Sec. 5]. 

A fundamental case in synchronization stability concerns with output-coupled 
identical linear systems under fixed interconnection. The problem is considered 
in [15| for time-invariant discrete-time systems and in [16] for continuous-time 
systems (as a generalization of Luenberger observer) leading to the following 
result: "If an individual system is detectable from its output and its system 
matrix is neutrally stable, then there exists a linear feedback law under which 
the trajectories of the coupled replicas of the individual system exponentially 
synchronize provided that the (directed) graph representing the interconnection 
is connected." We emphasize that (i) the result needs no assumption on the 
strength of coupling and (ii) synchronizing feedback law is independent of the 
number of systems and their interconnection. In this paper we extend this result 
for time- varying linear systems. 

For a time- varying pair (C, A), where A(-) is the system matrix and C(-) 
is the output matrix, we first define synchronizability (with respect to set of 
all connected interconnections.) Roughly, a pair (C, A) is synchronizable if one 
can find a bounded time- varying linear feedback law L(-) under which the tra- 
jectories of the coupled replicas of the individual system described by triple 
(C, A, L) synchronize for all connected interconnections. Then we study the 
conditions that would imply synchronizability. The assumptions and results 
almost parallel the time-invariant case. The assumption we make on the sys- 
tem matrix is that its state transition matrix is bounded in both forward and 
backward time0, which yields (considering trajectories) sustained and bounded 
oscillations. Boundedness in forward time is necessary for stability because we 
make no assumption on the strength of coupling. Boundedness in backward time 
can be relaxed at the expense of complicacy of analysis and need for additional 
technical assumptions on pair (C, A). For simplicity, therefore, we choose to 
keep it. One of the contributions of this work are in establishing the following 
results: 

• If pair (C, A) is asymptotically observable then it is synchronizable. 

1 Bounded state transition matrix assumption can be encountered in seemingly different 
problems in the literature whenever observability is at stake; see, for instance, [5]. 
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• If pair (C, A) is uniformly observable then it is exponentially synchroniz- 
able| 

Asymptotic observability we define as that the integrand of the observability 
grammian satisfies a general (yet technical) condition. This condition, which 
we name sufficiency of excitation H, is significantly weaker than persistence of 
excitation and allows the following result, cf. [T21 Thm. 2.5.1]. 

• Let Q be bounded and Q(t) — Q(t) T > for all t > 0. Linear system 
x = —Q(t)x satisfies lim^oo x(t) = if Q is sufficiently exciting. 

Uniform observability, on the other hand, is quite a standard concept, which is 
more or less equivalent to that the integrand (summand) of the observability 
grammian is persistently exciting. 

To obtain the above listed results we first study synchronization stability of 
a particular type of array. This array is pertinent to consensus problem (for 
trajectories are static once synchronized) yet different from the usual array of 
interest in consensus problems [3]. Our second contribution in this paper is in 
analyzing this new type of consensus array and, consequently, unraveling two 
arrays' similarities and differences. In addition, we investigate the stability of 
their union. The array dynamics generally studied in consensus problems is 

p 

xi = ^2nfij(t)(xj -Xi) (i) 

where Xi £ M. n is the state of the ith system and "fij(t) > for all t. What is 
known about this array is that its trajectories are bounded. In fact, the convex 
hull of the states co{xi, . . . , x p } is forward invariant regardless of the evolution 
of 7y (-)- Moreover, if certain connectedness property is satisfied by the graph 
described by {7^}, then trajectories Xj(-) meet at some common point, i.e. 
reach consensus. Finally, in general, there does not exist a quadratic Lyapunov 
function to establish stability; so the convex hull of the states is used instead 
[9] . The array considered in this paper is 

p 

Xi = ^2 ^viVi - Vi) > Vi = Q(t)%i ( 2 ) 

where time- varying output matrix Q(-) is symmetric positive semi-definite and 
7y is fixed. Below we list our findings residing in Section [4] 

• Like array ([1]), trajectories of array ([2]) are bounded. 

• Unlike array ([1]) , there exists a quadratic Lyapunov function^ for array ([2]) . 

2 Along with these results, we also provide a synchronizing feedback law L(-) in the paper. 
3 See Definition [2] 

4 However, the convex hull is no longer forward invariant. 
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• For Q sufficiently exciting, trajectories of array @ reach consensus for all 
connected interconnections. The point of consensus is independent of the 
evolution of Q. 

We also look at the union of the two cases Xi — Ysj=i 7y (*) (Vj ~ Vi)> Vi = 
Q(t)xi. We find that unbounded trajectories may result from this situation, 
hence stability is no longer guaranteed. 

The outline of the paper is as follows. After introducing notation and basic 
definitions, we define synchronizability and give the formal problem statement 
for continuous-time linear time- varying systems in Section[3J This section also is 
where we draw the simple link between synchronization of time-varying linear 
systems and consensus of array © . In Section 0] we establish the stability 
of array @ via a quadratic Lyapunov function and construct (observability) 
conditions on Q yielding consensus. In Section [5] we interpret these conditions 
through the observability grammian of time-varying pair (C, A) and establish 
our main results. Finally, in Section [SJ we generate the discrete-time versions 
of the continuous-time results. 



2 Notation and definitions 

Let N denote the set of nonnegative integers and R>o the set of nonnegative 
real numbers. The meaning of N>k is the obvious. Let | • | denote (induced) 
2-norm. Identity matrix in R" xn is denoted by /„. The set of all symmetric 
positive semi-definite (SPSD) matrices in R™ xn is denoted by Q n . We also 
define Q n := {R G Q n ■ \R\ < 1}- Let 1 G W denote the vector with all 
entries equal to one. The smallest and largest singular values of A G R mxn are, 
respectively, denoted by a m i n (A) and a max (A). Kronecker product of A G M. mxn 
and B G W xq is 



A®B 



a\\B ■ ■ ■ a ln B 



Q*rn\B • • • (l mn B 



Kronecker product comes with the following properties: (A <g) B)(C ® D) = 
{AC)®(BD) (provided that products AC and BD are allowed); A®B+A®C = 
A® {B + C) (for B and C that are of the same size); and (A®B) T = A T ®B T . 
Moreover, the singular values of (Ai&B) equal the (pairwise) product of singular 
values of A and B. 

A (directed) graph is a pair (Af, £) where Af is a nonempty finite set (of nodes) 
and £ is a finite collection of ordered pairs (edges) (rij, nj) with n,;, rij G Af. A 
directed path from n\ to ni is a sequence of nodes (m, ri2, ■ ■ • , rig) such that 
(m, rii+i) is an edge for i G {1, 2, ...,£— 1}. A graph is connected if it has a 
node to which there exists a directed path from every other node0 The graph 



5 Note that this definition of connectedness for directed graphs is weaker than strong con- 
nectivity and stronger than weak connectivity. 
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of a matrix M := [my] G R pxp is the pair (Af, £), where AT = {ni, ri2, . . . , n p } 
and £ is such that (n i} nj) £ £ iff my > 0. Matrix M is said to be connected 
when its graph is connected. 

Throughout the paper T :— [7y] G M pxp will represent an interconnection 
(in the continuous-time sense) satisfying 7y > for i 7^ j and 7^ = — Sj^i 7«j 
for all i. It immediately follows that A = is an eigenvalue with eigenvector 1, 
that is, n = 0. For V connected, eigenvalue A = is distinct and all the other 
eigenvalues have real parts strictly negative. Let rgP satisfy 

r T r = (3a) 
r T l = 1. (3b) 

Then r is unique (for T connected) and satisfies lim^oo e rt — lr T . Also, r has 
no negative entry. We denote by G>o the set of all connected interconnections, 
i-c G>q — {T G R pxp : T connected interconnection, p = 2, 3, . . .}. 

Matrix A := [Ay] G R pxp denotes an interconnection (in the discrete-time 
sense) satisfying Ay > for all i, j and J2j = 1 f° r au It follows that 
A = 1 is an eigenvalue with eigenvector 1, that is, Al = 1. For a connected A, 
eigenvalue A = 1 is distinct and all the other eigenvalues lie strictly within the 
unit circle. Let r£l p satisfy 

r T A = r T (4a) 
r T l = 1. (4b) 

Then r is unique (for A connected) and satisfies lim^oc A fe = lr T . Also, r 
has no negative entry. By slight abuse of notation (yet with a negligible risk of 
ambiguity) we will let Q >n also denote the set of all (discrete-time) connected 
interconnections A. 

Let § G {M> , N}. Given maps & : S -» R n for i = 1, 2, . . . , p and a 
map £ : § — > R™, the elements of the set {&(•) : i = 1, 2, . . . , p} are said to 
synchronize to £(•) if |&(s) — £(s)| — > as s — > 00 for all i. They are said to 
synchronize if they synchronize to some £(■). Moreover, if there exists a pair of 
positive real numbers (c, a) such that max^ |&(s) — £(s)| < ce~ QS for all s, then 
&(•) are said to exponentially synchronize. 



3 Problem statement 

For a given interconnection T = [7y] G R pxp , let an array of p linear systems 
be 

±i = A(t)xi + m (5a) 
y t = C(t) Xi (5b) 

Zi = 'Y^liAVj-yi) (5c) 

where Xi G R" is the state, Ui G M™ is the input, yt G R m is the output, and 
Zi G R m is the coupling of the ith system for i = 1, 2, . . . , p. For each t G R we 
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have A(t) € R nxn and C(t) £ R mxn . The solution of ith system at time t > 
is denoted by Xi{t). We denote by $a( - > •) the state transition matrix for A, i.e. 
the unique solution of the matrix differential equation 

with $^4(io, io) = In- Also, recall that the observability grammian of pair 
(C, A) is given by 

W o {t , *):=/ io)C T (r)C(r)^(T, t )dr (6) 

•/to 

for to, i £ t. We will henceforth assume that the integrand of the grammian is 
Riemann-integrable. 

Definition 1 (Synchronizability) Given functions A : R — » R" x " and C : 
M — >• R mx "; pair (C, A) is said to be synchronizable (with respect to G>o) if 
there exists a bounded, time-varying linear feedback law L : R — > R™ xm s«c/i 
that for each V G ^ >0 , solutions aJj(-) o/ array ([5]) wif/i it, = L(t)zi synchronize 
for all initial conditions. 

Our objective in this paper is to find sufficient conditions on pair (C, A), 
in particular on the observability grammian ([6]), for synchronizability and to 
design a synchronizing feedback law L when proposed conditions are met. 

The above statement of our objective almost suggests that we first find 
sufficient conditions and search for an L only afterwards. However we adopt 
the opposite approach. We choose first to construct an L and then work out 
the conditions on (C, A) for synchronization under such feedback law. Given 
(C, A) let 

L(t):=$ A (t,0)$>l(t, 0)C T (t). (7) 
For interconnection V £ R pxp consider array (O with ui = L(t)zi- We can write 
ii = A(t)xi + L(t)C{t) Hi Oj (8) 

Let us define the auxiliary variable £ 1" as 

Ut) := $ A (0, t) Xi (t) (9) 
for % = 1, 2, p and t > 0. Combining (J7J), ©, and © we obtain 

6 = 0)C T (t)C(t)$ A (t, 0)£ 7 y& ~ 6) • (10) 

Now note that if <& A is bounded, then synchronization of solutions £,;(•) implies 
synchronization of solutions Xi(-) by ([9]). Moreover, if C is bounded as well, 
then boundedness of L is guaranteed by ([7]). Based on this simple observation 
let us write the following assumption to be invoked later. 
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Assumption 1 (Boundedness) For A : K -> R" xn and C : M -> M mx " 
following hold. 

(a) There exists a > 1 suc/i £/ia£ |$^(fi, ^2 ) | < 5 /or a/7 ti, £2 > 0. 
(7>J There exists c > 1 suc/i t/ia£ |C(t)| < c /or all t > 0. 

Remark 1 Note that in the time-invariant case Assumption \^b) comes for 
free and Assumption^ a) boils down to that matrix A is neutrally stable (in the 
continuous-time sense) with all its eigenvalues residing on the imaginary axis. 

The second point we want to make is that the term multiplying the sum in 
(fT0|) is the integrand of the observability grammian, which is SPSD at each t. 
We elaborate on this fact in the next section. 



4 Synchronization under SPSD matrix 



For a given interconnection T = pyy] 6 R pxp , let an array of p systems be 



E 



7y j 



(11) 



where %i € M" is the state of the ith system (for i = 1, 2, . . . , p) and Qt € ffi nxrl 
is SPSD for each t > 0. We assume Q : K>o - * Qn to be Riemann-integrable. 
By letting 



x := 



we can rewrite more compactly as 

x = (r®Q t )x. 



(12) 



Remark 2 Sometimes we need function Q : R>o — > Q n 6e bounded on the 
interval [0, oo), i.e. f/iere exists h > 1 si/c/i i/iai |Q t | < ft, /or aZZ iVoie i/ia£ 
(fT5|) can 6e written as 



Now, since T is an interconnection, so is hT. Also, connectedness is invariant 
under multiplication by a positive scalar, i.e. T is connected if and only if hT is. 
Finally, observe that Qt/h G Q n - Without loss of generality (for our purposes) 
therefore we can can take h to be unity, which lets us consider Q : R>o — » Q n 
whenever we need Q be bounded. 
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In the rest of this section we first show that the origin of system (fT"2"]) is 
stable regardless of interconnection T £ M. pxp and function Q : R — > Q n . Then, 
under connectedness of T, which is obviously necessary for synchronization, we 
work out some sufficient conditions on function Q to establish synchronization 
of solutions Xi(-) of array (fTTj) . Finally, we provide two theorems to make the 
picture that we want to give in this section closer to complete. One of those 
theorems states that time-invariance of interconnection T in (|12|) is necessary 
for stability. With the other one, we aim to show that the sufficient conditions 
that we will have proposed on Q for synchronization cannot be readily relaxed 
into a less technical one without sacrificing generality. 

4.1 Stability 

Lemma 1 Let interconnection T £ W xp be connected and r £ R p satisfy (|3j). 
Then, there exists symmetric positive definite matrix £1 £ W xp such that 

(r - lr T ) T n + fl(T - lr T ) = -I p . (13) 

Proof. Consider matrix T — lr T . Observe that (r — lr T ) k = T k + (-l) fc lr T 
for k £ N. For igK, therefore we can write 

e( r - lrT )* = I p + t(T-lr T ) + j(T-lr T ) 2 + ... 

I p + tT + l -Y 2 + . . )j - (tlr T - f -lr T + . 

= e rt -(l-e^)lr T . 

Consequently, lim^oo e^- r ~ lr = 0; and we deduce that [r — lr T ) is Hurwitz. 
Therefore Lyapunov equation p"3|) admits a symmetric positive definite solution 

n. m 

Lemma 2 Let interconnection T £ M. pxp be connected, r £M. P satisfy |3]), and 
symmetric positive definite matrix Q £ W xp satisfy (| 13|1 - Define V : R" p — > 
R> as V(x) := x T (fi (g> 7 n )x. Then, for all Q : R> -> Q„ and all t > 0, 
solution of system (I12p satisfies 

jV(x(t) - x) = -(x(t) - x) T (/ p ® Q t )(x(t) - x) 

where x := (lr T ® J n )x(0). 

Proof. Observe that (lr T ® J„)x(t) = 0, which implies (lr T ® J„)x(t) = x for 
all i > 0. Whence x(t) = ((T - lr T ) (g> Q t )( x W ~ x )- We can therefore write 

jV(x(t)-5L) = x(i) T (fi®/„)(x(t)-x) + (x(t)-x) T (r!®/„)x(t) 

= (x(i) - x) T ((r - lr T ) ® Q t ) T (ft ® in)(x(i) - x) 
+(x(t) - 5c) T {n ® J n )((r - lr T ) ® Q t )(x(i) - x) 
= (x(i) - x) T [((T - lr T ) T VL + 0(r - lr T )) ® Q t ] (x(i) - x) 

= -(x(t)-x) T (/ p ®g t )(x(t)-x). 



8 



Hence the result. 



Theorem 1 (Stability) Given interconnection T £ R pxp , there exists a > 
such that, for all Q : R>o — > 2n, solution of system (JT3J) satisfies 



|x(t)| < a|x(0)| 



/or a// 1 > 0. 



Proof. Interconnection T is similar to a block diagonal matrix diag(ri, F2, . . . , T q ) 
in M. pxp such that Ti £ M. PiXPi for i = 1, 2, . . . , q is a connected interconnection 
if Pi > 2 and I\ = otherwise. (Integer g equals the number of eigenvalues of T 
at the origin.) Since diag(r i7 T 2 , . . . , T q )®Q t = diag(ri<g)Q t , r 2 <g>Qt, ■ ■ ■ , L g <g) 
Qt) without loss of generality it suffices to check two cases: (i) T = 0; and (ii) 
T is connected. First case is trivial; so let us suppose F is connected. 

Now, let r £ M. p and symmetric positive definite matrix fl £ R pxp satisfy 
([3]) and (TT3]) . respectively. Given Q : M>o — > Q n , consider system (fl2|) . Let 
x = (lr T eg) J n )x(0). Recalling that r has no negative entry we can write 

|x| < |lr T ||x(0)| 
< |l||x(0)| 
= V^|x(0)|. 

Let V(x) = x T (f7 (g) i„)x. Lemma [5] yields that ^(x(-) — x) is nonincreasing. 
Hence 

|x(i)| < |x(i)-x| + |x| 

< . 1 , n y nx(t)-x) + |x| 

1 



< . — _ Vy(x(o)-x) + |x| 

V 0"min(") 

< /£^| x(0 )-x| + |x| 



< 



^^a (1+VS ) + vs)ix(o)i 



for all t > 0. ■ 

Theorem[T]establishes stability. That is, for a fixed interconnection^] T, which 
need not be connected, solutions Xi(-) of array Ijlip stay within a bounded region 
(that depends only on initial conditions xi (0) and interconnection T) for all time- 
varying SPSD matrix Q. Whenever T is connected, that bounded region can be 
described quite precisely. See the next result, which is a direct consequence of 
Lemma [2] 

6 Later in the section, we will also investigate whether stability is preserved when both T 
and Q are time-varying. 
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Theorem 2 Given connected interconnection r £ K pxp , let r 6 MP and sym- 
metric positive definite matrix Q, £ M. pxp satisfy Q and (fT3|) , respectively. 
Then, for all Q : M>q — > Q n , solutions Xi(-) of array (jlip satisfy, for all t > 0, 



where x := (r T ® 7„)x(0). 

4.2 Asymptotic synchronization 

We now begin looking for sufficient conditions on Q : R — > Q n that guarantee 
that solutions ajj(-) of array (jTTJ) synchronize. The next fact is to be used by 
the key theorem following it. 

Fact 1 Let f : [0, T] — > [0, 1] fee Riemann-integrable. Then 



Proof. Result trivially follows for T = 0. Suppose T > 0. Fix some <5 > 
such that T/6 =: N is an integer and let If := Jj Q „ /. Then there exists 
fci 6 {1, 2, ... , TV} and tx £ [(fci - 1)<5, faS] such that /(ii) > I f /T. Since 
/(i) _■ 1 f° r a h i, we can also claim that there exists fc 2 € {1, 2, ... , iV} \ 
{fci} and t 2 G [(fc 2 - fc 2 <5] such that /(i 2 ) > (If - 8)/T. Following the 
pattern, we can generate a sequence (ki)\!j-/ S of distinct elements from the set 

{1,2, . . . , N} and an associated sequence satisfying t{ £ [(hi — 1)5, kiS) 

and f(ti) > (// — (i — 1)6) /T. Therefore we can write 




1/2 




LT/aj 



E sup J 2 (i)<5 > 



fc=1 te[(fc-l)5, fci] 
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By definition of integral, we can therefore write 

,T l T / s i 

/ f 2 (t)dt = lim V sup f{t)8 

JO s ^0+ • ' tGUk-l)8.k8] 



I f -(i- 1)5 



6^0+ ^ V T 

1=1 



k=l te[(fe-l)5, kS] 

llf/i 

> lim S V 

! = 1 
r U>/«J 

2= 1 

- A ™ E (3 - 2I f* s ) + E ^ 



i=l i=l 



= lim — =r « 2 

i=l 

Hence the result. ■ 
Theorem 3 Given a pair of positive real numbers (e, T), define 

S(e, T) := mini-, & A . (14) 

Given connected interconnection T S R pxp , Zet r € R p and symmetric positive 
definite matrix Q € R pxp satisfy ([3]) and (|13p . respectively. Define 

p(r) := f max 

(fi)max{l,|r| 3 }. (15) 

Let 7(x) := x T (fi <8> 7„)x /or x e M™p. T/ien, /or aZ/ Q : M> -> Q„ } the below 
inequality 

O-min ^jf Qtdtj > e (16) 
implies that solution x(-) of system (|12p satisfies 

V(x(T) - x) < (l - ^Jp) K(x(0) - x) 
where x := (lr T ® J n )x(0). 

Proof. Given pair (e, T) let w := e/(4T). Consider system (|12l) . Let us 
introduce 

:=x(i)-x. 



L^//5J 
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By Lemma [2J we have 

vm) = -t T (w P ®Qt)m- 

Also, £(•) can be shown to satisfy 

£ = (r®Q t )£. 



(17) 
(18) 



Let Q : R>o — > Q„ satisfy (116|) . Then, regarding the evolution of £(•), one of 
the two following cases must be. 

Case 1: \((t) - £(0)| < w|f(0)| for all i G [0, T]. Let &(i) := - £(0) 
and recall that |Qt| < 1. For reasons of economy, let us adopt the notation 
Qt := (I p ®Qt)- Note that then we have |Q t | < 1 as well as 



Qtdtj > e . 

(Observe that T > e.) From ([T7]) we can write 



vim) = mm 



£ T (t)QtS(t)dt 







[ Q t dt 


£(0)-2 / 


Jo 


Jo 



< ne(o))-e T (o) 



< vm)-e\m\ 2 + ^T\ao)\ 2 



b T (t)Q t b(t)dt 



< i 



2cr max (n) 



(19) 



Case 2; |£(*)-f(0)| = w|£(0)| for some i £ (0, T]. Without loss of generality, 
assume |f (i) - f (0)| < w|£(0)| for t e [0, F). We can by |Q2J) write 

f\Qti{t)\dt = irr 1 / |r<8J n ||Q t e(t)|tft 

Jo Jo 



> irr 1 

> irr 1 



l(r®Q t )$(t)|* 



(r Q t )f (t)A 



|rr^|e(o)| 



(20) 
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Since |Q t £(t)| < (l+w)|£(0)| for t € [0, t\, we can invoke Fact Q] on (2© and 
obtain 



> (1 +1 ,) 2 |«0)| 3 



\ 3 1 

W \ 1 



.(l + wjiri; 3P 

> (i + w ) 2 |c(o)| 2 ' 



1 



3(l+w)T 2 |r| 3 



(i + w)|r|y 3r 2 

|£(0)| 2 . (21) 



1 /2 

Since |Q t | < 1, we have Q/ > Q t . Also, note that V(£(t)) is nonincreasing 
thanks to (jTTJ) . Now, by (|2"Tj) we can write 



T 



V(£(T)) = V(f(0)) - / e J (t)Q t £(t)dt 

Jo 

Jo 

= vm)- f n\ ,2 m\ 2 dt 

Jo 

< V(f(0))- [ \Q t m\ 2 dt 



,3 



* ^ »- 3(i + ,)T 2 |r|3 ^)l 2 



„3 



5 ( 1 -w^b«)) , '~- (22) 

The result follows by ([13 ]) . ([15 ]) . ([13 ]) , and ^2 ]) . ■ 

Fact 2 Lef {a,i)^L 1 be a sequence with < a% < 1 /or aZ/ i. Then, product 
— a converges to zero if and only if sum a i diverges. 

Proof. See e.g. [12j Thm. 1.17 of Ch. VII]. ■ 
In the light of Theorem[3]and Fact[2J we now state our most general condition 
on Q : R — > Q n for synchronization of solutions Xi(-) of array (fTTj) . 

Definition 2 Function Q : R>o — » Qn *s to fre sufheiently exciting if there 
exists a sequence of pairs of positive real numbers (e,-, Ti)^. 1 satisfying 

ti+Ti \ 

Qtdt > e< (23) 



13 



for ti = YTj—i Tj with t\ = 0, and 



N 

N—too 



lim V<5( £l , Ti) = oo (24) 

/ — >nn ' ■ 



where <5(-, •) is as defined in (| 14|) . 

Theorem 4 Let Q : R>o — > Q n &e sufficiently exciting. Then, for all connected 
interconnection T S M pxp ; solutions Xi(-) of array (fTTj) synchronize to 

x{t) = (r T ®/„)x(0) 

where reP satisfies ([3]). 

Proof. Let us be given function Q : ^>o - ► Q n an d sequence (e^, Ti) satisfying 
P^j) and (El- Let us let ^ := 5{e.^T l ). Let interconnection T e M pxp be 
connected, r € R p satisfy ([3]), and symmetric positive definite matrix f2 € R pxp 
satisfy (fT3|) . Consider system (fT2|) and let x = (lr T ® 7„)x(0). Now, letting 
V^x) = x T (f2 £g> Jrt)x, by Theorem [3] we can write 



v(x(t i+ i) - x) < \i - -^j y(x(t0 - x) 

for i = 1, 2, ... where p(-) is as defined in (fTS)) . Whence 



i-l 



7(x(ti) - x) < V(x(0) - x) J] ( 1 - J ■ (25) 



3=1 

Let now a.; :— Si/p(T). We can write by (j2~4")) 

iV iV 



lim >^ a.; = p(r) 1 lim >^ 5i 

i=l i=l 



Now, we can invoke Fact [5] and claim that limjv^oo ]X =1 (1 — a.;) = 0, which 
yields by 



lim F(x(t;)-x) = 0. 

i— >oc 

Hence the result. ■ 

Corollary 1 Let Q : R>o — > Q n &e sufficiently exciting. Then solution to linear 
system x = —Qtx satisfies lim^oo a;(f) = 0. 

The following definition is quite standard; see, for instance, j8j. Note that 
the condition it depicts is less general than that of Definition^ yet it guarantees 
exponential synchronization. 
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Definition 3 Map Q : M>o — > Q n is said to be persistently exciting if there 
exists a pair of positive real numbers (e, T) such that 



t+T 



rdr > e 



(26) 



/or aZZ t > 0. 



Remark 3 The following theorem can be viewed as a generalization of a classic 
result in adaptive control theory USH Thru. 2.5.1]. Note that Corollary]]] makes 
another generalization to this result since "sufficiently exciting" is weaker than 
"persistently exciting". 

Theorem 5 Let interconnection T G M. pxp be connected and function Q : 
K>o — * Q n be persistently exciting. Then, solutions Xi(-) of array ex- 
ponentially synchronize to 

x(t) = (r T ®/„)x(0) 

where r £ K p satisfies ([3]). 

Proof. Since Q : R>o — > Qn is persistently exciting, by definition, there exists 
a pair of positive real numbers (e, T) satisfying (|26|) for all t > 0. Let V(x) = 
x T (f2 (g) J n )x for x S M np , where symmetric positive definite matrix f2 £ R pxp 
satisfy (| 13[) . Then, by Theorem [3] we can write 



V(x(fcT) - x) < 1 



p(r) 



K(x(0) - x) 



for all k £ N, where x := (lr T ® J„)x(0). The result then follows. ■ 
We now present an interesting application of Theorem[5]on coupled harmonic 
oscillators (in M 2 ) described by 



Xil 
Xil 



Xt2 



-xn + ^2lij{Xj2 - x i2 ) 



for i = 1, 2, 



p. Let 



Xi2 



A 



1 
-1 



, C=[0 1] 



Then we can write 



Xi = Aa; 4 + C T C ^ 7y (xj - x{) . 

Define &(t) := e - At ii(t) and £ = [£f ... £j] T . Then 

C= (r®e ATt C T Ce At K 
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for A is skew-symmetric. A trivial computation shows 



a 1 twr 



C 1 Ce 



A I 



— sin t cos t 



sin t cos t 
sin 2 i 



6 Q 2 



whence 



t + 27T 



e A T r C T Ce Ar dr 



7T/ 2 



for all t. Therefore t ^ e A t C T Ce At is persistently exciting. Now, suppose 
that r is connected and r S M p satisfies ((3|). Then, by Theorem [5] solutions 
exponentially synchronize to = (r T ® I 2 )£(0). Since Xj(t) = e At £i(t) and 
e' 4 * is an orthogonal (hence norm-preserving) matrix, it follows that solutions 
Xi(-) of the coupled harmonic oscillators exponentially synchronize to 

"zi(Q)' 



x(t) = (r 1 



a?p(0) 



4.3 Negative results 

Before we end this section, we present two negative results, which we believe 
constitute answers to questions that arise naturally. The first of those questions 
emerges as follows. In Theorem [T] we have proven that system x = (T eg) Qt)x, 
where T is a fixed interconnection and Q is a time-varying SPSD matrix, has 
a bounded solution for all initial conditions. It also trivially follows from the 
results in, for instance, O [7] that solution of system x = (r t ® Q)x, where this 
time interconnection T is time- varying and SPSD matrix Q is fixed, is bounded. 
At this point, it is tempting to ask the following question. 

Is solution of system x = (T t <g) Qt)x bounded? 

The answer is not always and it is formalized in the below result. 

Theorem 6 There exist maps T : M — ► M px?) , where T t is an interconnection 
for each t, and Q : R — ► Q n such that system x = (r t Qt)x has an unbounded 
solution. 

Proof. We construct T and Q as follows. Let 

" 

0- 1 1 

1 0-10 


" 

1- 1 
0-11 




r a := 



-1 


1 






































1 


-1 


-1 





1 
































1 





-1 



1G 



and 



" 1 


" 




'00" 




' 0.5 


0.5 " 


, Qd := 








, Qb ■ — 


1 


, Qc := 


0.5 


0.5 



Now let both T and Q be periodic with period T = 40 with 

Q 



0.5 
-0.5 



r, 



r a for < t < 10 

T b for 10 < t < 20 

T c for 20 < t < 30 

T d for 30 < t < 40 



and Q t 



for < t < 10 

Q b for 10 < t < 20 

Q c for 20 < t < 30 

Q d for 30 < t < 40 



Whence we can write for k — 0, 1, 



x(fcT) = A fc x(0) 



where A := e (r d ®Qd)io e (r c ®Q c )io e (r b ®Q b )io e (r a ®Q a )io^ when the eigenvalues of 
A are numerically checked, one finds that there is an eigenvalue outside the unit 
circle (|A| w 2), which lets us deduce that the origin of system x = (T t ® Qt)x 
is unstable. I 
Our first question was concerned with stability under time-varying inter- 
connection; and we have seen that solutions Xi(-) of array (| 1 1 j) need not stay 
bounded in such a case. The second question is about synchronization. By 
Theorem 2] we know that if map Q : IR>o — * Q n is sufficiently exciting, then for 
all connected interconnection T solutions of array Ijlip synchronize. Now, sup- 
pose that we are given some sufficiently exciting Q with an associated sequence 
( £ ii Ti)°°^ ll see Definition [2] Note that, due to (fT4|) . we have £ i = °°- I n 

addition, since Qt € Q n for all t, we have Ti> Ei, which yields Y^Li Ti — oo. 
Applying these observations on ([2^)1 . we obtain 



-0.5 
0.5 



lim a n 



Qtdt 



oo . 



(27) 



Above condition, depicted in (|2T|) . can be shown to be necessary for synchro- 
nization. Now, we ask the following question. 

Is condition (j27[) sufficient for synchronization ? 

The answer turns out to be negative. In fact, even a much stronger condition 
is not sufficient for synchronization as the following result shows. 

Theorem 7 There exist connected interconnection T 6 M. pxp and map Q : K — > 
Q n satisfying 



1 



lim inf — cr min 

T — *oo 1 




ttdtj > 

such that solutions Xi(-) of array (|11| do not synchronize. 
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Proof. As in the proof of the previous result, we will once again make use of 
projection matrices. Let Ek = 2~ fc for k = 1, 2, . . . Then, let piecewise linear 
function <d : R> -> R be 

:= i9(r fc _i) +tsin£ fe cos£ fe for < e [r fe _x, r fe ) 

where r = 0, i9(0) = 0, and 

7"fe = Tk-l 



Sin £fc COS £fc 

0(t*) = lim - e k+1 



for fc = 1, 2, . . . Now define Q : M>o — + Q2 as 

cos 2 i?(t) smi?(t)costf(t) 



sint?(t) costf(t) 



sin 2 



which is a projection matrix that projects onto the line spanned by the vector 
[cos#(t) sin$(i)] T . Calculations yield 



1~k ~ Tk-l 



Qtdt = -h ■ 



Observe that 



and 



lim 



Tfc-l 1 



ah 



k — >oo Tk 



Qtdt = a + a ^ 



Qtdt 



for all a > and ti>t\> 0. Given T G (Tfe_i, Tfc] we can write 



Q t di = 



0"n 



> 
> 

which yields by flUJ that 

lim inf — er 1T 



T 
1 

T 

Tfc-1 

2T 



Tfc-1 



Qtdt 



Qt* > - . 



(28) 
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Let us now consider pi[) under the following connected interconnection 



T := 



-1 1 




Setting £2(0) = we can write 



(29) 



Note that we need limt_>oo x\ (t) — for synchronization since a^O) = 0. In 
terms of polar coordinates, i.e. x% = [rcos6 rsin6] T , we can express (|29[) as 



= -rsin 2 (tf(*) -9 -tt/2) 

= sin(i?(i) -6- tt/2) cos(tf(t) - - tt/2) . 



(30a) 
(30b) 



Let us initialize cci such that r > and 0(0) = —tt/2 — ej. We then observe 
that 6{t) = &(t) for all t > 0. Eq. ([3D]) simplifies to 

• 2 

r = — rsm 

= sin£fcCOS£fc 

for t € [Tfe—i, tj.) and fc = 1, 2, ... Thence 

r(r fc ) = r(r fe _ 1 )e-^- T *- 1 ) sin2£ " 



which yields 



"(Tk) = r(0) exp ( -2-7T ^ tan 

/ fe N 

= r(0)exp ( -27r^tan2~* 



> r(0) exp I 

> r(0)e- 47r 
for all fc. Therefore limt— >oo #1 (i) 7^ 0. 



5 Observability grammian and synchronizabil- 
ity 

Based on the results of the previous section, we are now ready to establish our 
theorems that are aimed to reveal the correlation between synchronizability and 
observability grammian. We begin with two definitions. 
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Definition 4 For A : R -> R nx ™ and C : R -> R mxn ,. pair (C, A) is said 
to fee asymptotically observable if the integrand of the observability grammian, 
1 1-> 0)C T (t)C(t)$ J 4(i, 0), is sufficiently exciting. 

The following definition is borrowed (with slight modification) from [6]. 

Definition 5 For A : R -> R nx ™ and C : R -> R mxn , pair (C, A) is said to fee 
uniformly observable if i/iere exists a pair of positive real numbers (e, T) smc/i 
i/iai 



^"min 

/or a// i > 0. 



Remark 4 For A : R -> R" xn and C : R -> R mxn satisfying Assumption^ 
asymptotic observability of pair (C, A) implies uniform observability of (C, A). 
For a time-invariant pair, which need not satisfy Assumption [7J' asymptotic 
observability, uniform observability, and the standard definition of observability 
(for time- invariant linear systems) are all equivalent. 

The next result is our main theorem. It states that a time-varying pair 
(C, A) is synchronizable if it is asymptotically observable. 

Theorem 8 Let A : R -> R nxn and C : R -> R mxn satojt/ Assumption^ If 
pair (C, A) is asymptotically observable, then it is synchronizable. In particular, 
if we choose L : R — » R nxn as in {7)1. i/ien /or eac/i T £ <?>o> solutions Xi(-) of 
array {5| u>i£/i u, = L(t)zi synchronize to 



x(t) := (r T ® $ A (t, 0)) 



si(0) 

Zp(0) 



(31) 



where r € R p satisfies ([3]). 

Proof. Let L be as in ([Jj; then it is bounded by Assumption Q] Let T € R pxp 
be a connected interconnection and r GMP satisfy ©. Consider array ([5]) with 
u.j = L(t)zi. Define auxiliary variables as 



Si(t)=* A (0, t)xi(t). 



(32) 



Then, we can write 

where Q t := (ac) _1 $^(t, 0)C T (t)C(t) < i> J 4(t, 0), 7y := acyy, and a, c > 1 come 
from Assumption [TJ Note that Q t & Q n for all t > 0, t i— > Q t is sufficiently 



20 



exciting, and T := [7^] is a connected interconnection satisfying r T T — 0. We 
now invoke Theorem [4] on p3|) to deduce that solutions synchronize to 



£(t) = (r T ®/„) 



i P (0) 



(34) 



Recall that $>a is bounded. Hence, combining (|34p and ([3^)1 yields (j3"Tj) . ■ 
As noted in Remark 31 uniform observability is more restrictive a condition 
than asymptotic observability. However, it has a stronger outcome as stated by 
the following theorem. 

Theorem 9 Let A : K -> R" xn and C : R -> R mxn satisfy Assumption^ 
Then pair (C, A) is synchronizable if it is uniformly observable. In particular, 
if we choose L : R — > R" xn as m 0, £/ien /or eac/i T G <?>o, solutions Xi(-) of 
array J5| mi/i u, = L(t)zi exponentially synchronize to 



x(t) := (r T ®$ A (t, 0)) 



ari(0) 
zp(0) 



where r € W satisfies ^ . 



Proof. The demonstration flows in a way that is analogous to that of Theo- 
rem [5] This time, however, the result follows from Theorem [5] ■ 



6 Discrete-time results 

Our study on synchronization has hitherto been solely in continuous time. How- 
ever, it is possible to extend the analysis to systems in discrete time without 
much difficulty. In fact, most of the continuous-time results have discrete-time 
counterparts under assumptions and definitions that are analogous to the ones 
we used for continuous-time systems. In this section, therefore, we focus on 
discrete-time time-varying linear systems and investigate the correlation be- 
tween synchronizability and observability grammian in discrete time. 

For a given interconnection A = [Ay ] € R pxp , consider the array of p discrete- 
time linear systems, for k G N, 

x~l = A{k)xi + Ui (35a) 
Vi = C(k) Xi (35b) 
Zi = ^K 3 {yj - Vi) (35c) 

where x% € R™ is the state, xf is the state at the next time instant, Ui £ R 71 is 
the input, y.- L £ R m is the output, and Zi € R m is the coupling of the ith system 
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for i = 1, 2, . . . , p. For each k £ N, we have A(k) £ R nxn and C(k) £ E. mxn . 
The solution of ith system at time k £ N is denoted by Xi(k). We denote by 
&a(', •) the state transition matrix for A, i.e. for k > ko 

$ A (fc, ko) = A(k - l)A(fc - 2) • ■ ■ A(fc ) 

with <f>yi(fco, feo) = In- We will let <&ji(ko, k) = fco) whenever the inverse 

exists. Observability grammian for pair (C, A) is given by 

fc-i 

W o (k , k) := ]T ^oJ^fflCf^.ff, k ) 

e=k 

for fc, fco £ N. Below we provide the discrete-time versions of Definition Q] and 
Assumption [TJ 

Definition 6 (Synchronizability) Given functions A : N — > M" x ™ and C : 
N — * R mx "; pair (C, A) is said to be synchronizable (with respect to G>o) if 
there exists a bounded, time-varying linear feedback law L : N — > M™ xm such that 
for each A £ <?>o, solutions Xi(-) of array (|35p with Ui = L{k)zi synchronize for 
all initial conditions. 

Assumption 2 (Boundedness) For A : N -> R nxn and C : N -> M mxrl 
following hold. 

(a) For each k £ N , A -1 (fc) exists. There exists a > 1 suc/i i/wrf &2)| < 

a /or aZZ fci, &2 € N. 

There exists c > 1 swc/i i/iai |C(fc)| < c /or aZZ £ N. 

Remark 5 When A and C are constant matrices, Assumption\^(b) comes for 
free; and Assumption\^(a) becomes equivalent to that all eigenvalues of matrix 
A are with unity magnitude and none of them belongs to a Jordan block with 
size two or greater. 



6.1 Synchronization under bounded SPSD matrix 

This subsection will emulate Section [U where we studied the stability and syn- 
chronization properties of array (jlip . For a given interconnection A = [Ay] £ 
W xp , let an array of p systems be 

xf = Xj + Qk Ay (asj — Xj) (36) 

where Xi £ W 1 and Qfc £ Q n for all fc £ N. We consider (|36)) as the discrete-time 
analogue of (fTT|) . Let us stack individual vectors into 
Then we obtain from (|36|) 

x+ = (7 np + (A - J„) <g> Q fc )x (37) 

which makes the analogue of system (fT2|) . 
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Remark 6 Recall that to establish stability of (continuous-time) array (jTTJ) it 
sufficed that Q t is SPSD for each t. (See Theorem^) That is, boundedness 
was not required. Later, when we established synchronization in Theorem [7J 
we needed solely that Q : K — > Q n is bounded. (See Remark^) The story 
has to be a little bit different in discrete-time. Note that in (|36[) we stipulated 
that Qk (the discrete-time counterpart of Qt) be in Q n . Even to be able to 
establish stability, let alone synchronization, we will need \Qk\ < 1, a more 
restrictive condition than boundedness. Clearly, this has to do with the fact that 
for discrete-time systems the magnitude of the righthand side is important for 
stability; whereas in continuous time, what matters (for stability) is only the 
direction of the righthand side. 

Lemma 3 Let interconnection A € W xp be connected and r € W satisfy (|4|). 
Then, there exists symmetric positive definite matrix £1 £ M. pxp such that 

(A - lr T ) T ft(A - lr T ) - O = -I p . (38) 

Proof. We first observe that (A — lr T ) k = A k — lr T . Then we can write 
linife^oo A fe — lr T = 0, which implies that matrix [A — lr T ] is Schur, i.e. all of 
its eigenvalues are strictly within unit circle. Therefore, discrete-time Lyapunov 
equation (|38|) admits a symmetric positive definite solution f2. ■ 



Lemma 4 Let interconnection A 6 W xp be connected, r S R p satisfy (l4|) . and 
symmetric positive definite matrix Q S W xp satisfy (J38J) . Define V : R np — > 
K>o as V(x) := x T (S}®I n )x. Then, for all Q : N — > Q n and all k G N, solution 
of system (|37[) satisfies 

V{x{k + 1) - x) - V(x(k) - x) < -(x(fc) - x) T (/j, <g> Q 2 k )(x{k) - x) (39) 

where x := (lr T ® 7„)x(0). 

Proof. Observe that (lr T (g) I n )x(k + 1) = x(fc) whence (lr r (g> I n )x{k) = x 
for all fc € N. Let f := x - x and A := A - lr T . Then we have = 
(I nj) + (A - I p ) <g> Qk)£,- We can write 

l/(C + )-^(0 = e{{ln P + {K~I p )®Qk) T {^®In){ln P + {K-I p )®Qk)-^®f}ji 

= C T (j(A Q - / p ) T r! + fi(A D - I p )) <g> (Q fe - Q2) + (A^fiA - fi) ® Qlj £ . 

Note that Q fe - Q| > since |Q fe | < 1 and that (A D - I p ) T fl + fl(A - I p ) < 
since A^57A Q — £1 < (this is almost immediate when we recall that the sublevel 
sets of quadratic Lyapunov functions are convex surfaces.) Therefore 

v(z + )-v(t) <-e(i P ®Qi)Z- 

Hence the result. ■ 
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Remark 7 Were Qk a projection matrix, then inequality (I39p could be replaced 
by the below equality 

V(x(k + 1) - x) - V(x(fc) - x) = -(x(fc) - x) T (/ p <g) Q fc )(x(fc) - x) . 

Theorem 10 Given interconnection A £ R pxp , t/iere exists a > suc/i that for 
all Q : N — > Q n , solution of system (|57|) satisfies 

|x(fc)| < a|x(0)| 

/or a// fceff. 

Proof. Demonstration uses Lemma |4] and flows similar to that of Theorem [TJ 

■ 

Fact 3 For map Q : N — > Q„, rea/ number e > 0, and integer N > 1 we have 
/N-l \ /N-l \ 



'mill J Qfe J > £ =^ C m i n I Q 2 

,fe=0 / Vfe=0 



e 2 

> 



iV 2 n 2 



Proof. Let us be given any neK" with v v = 1. Let er min (X}fe=o QkJ > £ > 
0. We can write 

v t Qqv + v T Qiv + . . . + v T Q N ^iv > e 

which implies that there exists k* £ {0, 1, . . . , N— 1} such that v T Qk*v > e/N. 
Since Qk* is an SPSD matrix there exists an orthogonal matrix R £ M. nxn and a 
diagonal matrix D £ R nxn with entries d, € [0, 1] for i = 1, 2, . . . , n such that 
Q fc , = R T DR. Also note that Q 2 t = R T D 2 R. Let [wi w 2 ■ ■ ■ u>n] T = w := 
Note that ui T w = 1 for i? is orthogonal. We can therefore write 

n 

which implies that there exists i* € {0, 1, . . . , n} such that di*w% > e/(Nn). 
Since u> 2 „ < 1, we can write 



£ 2 



A^ 2 n 2 



< 4<4 



< ^ £tf 

\k= 

whence the result follows, for v was arbitrary. 



< d%wl 

< w T D 2 w 
= v T Q\,v 

/N-l 

k 

fc=0 
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Theorem 11 Let e > be a real number and N > 1 an integer. Define 

Let interconnection A G M pxp be connected, r G K p and symmetric positive 
definite matrix €E IR pxp ; respectively, satisfy ((4| and ()38|) . Define 

p(A) := (ft)max{l, |A-/ P | 2 }. (41) 

Let F(x) := x T (fi (g) 7„)x /or x e M" p . Traen, for all Q : N -> Q„, tfie &eZow 
meguaZii?/ 

/JV-1 \ 



E Q* " 6 (42) 



implies that solution of system (I37[) satisfies 

V(x(N) - x) < (l - ^(x(0) - x) 

where x := (lr T (g) 7„)x(0). 

Proof. Given pair (e, T) let := e 2 /(AN 3 n 2 ). Consider system (1571) . Let us 
introduce 

£0) :=x(fc) -x. 

By Lemma HI we have 

V(£(k + 1)) - V(t(k)) < -f{k)(I p ® Q^)C(fc) . (43) 
Also, £ can be shown to satisfy 

e = (/«, + (A-Jp)®Q fc )^- ( 44 ) 

Let Q : N — > Q„ satisfy (|42| . Then, regarding the evolution of £(•), one of the 
two following cases must be. 

Case 1: |£(fc) - £(0)| < w|£(0)| for all fc G {1, 2, . . . , N - 1}. Let 6(fc) := 
£(fc) — £(0) and recall that \Qk\ < 1- Let Qfe := (7 P <g) Qk)- Note that then we 
have | Qfe | < 1 as well as, by FactO 

ffmin E Qfe = CTmin [ L p ® Q fe 

\fc=0 / \ k=0 J 

> ^ 

" N 2 n 2 ' 
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From (pt3|) we can write 

N-l 



fc=0 
N-l 

= v(t(o)) E (^(°) + Kfc)) T Q?-(C(o) + b(k)) 

k=0 

/N-l \ AT-1 AT-1 

= ^(e(o)) - e T (o) E q h ^ (0) - 2 £ feT ( fc )Q^(°) - E &T ( fc )Q^) 



> fe=0 / fc=0 fc=0 



^ ^(0))-^|£(0)| 2 + 2^(0)| 2 
= ^«°))-2^?l«°)l a 

=•2 



5 I'- iWMl ) 1 "- < 45 > 

Case 2: |f(jfc) - £(0)| > w|£(0)| for some jfe G {1, 2, . . . , N - 1}. We can by 
lgl|) write 



Eiq^wi = iA-/ P r i Ei(A-/ P )®/niiQ^(fc)i 

fc=0 fc=0 

k 

> lA-jpl^EKCA-^^Q^)! 

k=0 
k 

= \A-i p \-^\t(k+i)-ak)\ 



k=0 
k-1 



Ee(fc + i)-c(fc) 

fe=0 

= lA-jpi- 1 !^)-^)! 

> \A-I p \- x u\Z(0)\- (46) 

Eq. (|46|) implies that there exists fc € {0, 1, . . . , k} such that |Qfc£(fc)l > |A — 
I p \- l u}\£(0)\/N which implies 

N-l 

EiQ^(fc)i 2 >|A-/ P |-v|e(o)i 2 /iv 2 . 

fc=0 
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Then, by (|4"5|) we can write 



N-l 



V(Z(N)) = ^(0))-^^(fc)Q^(fc) 

k=0 

< ^(°))- F ^ie(0)l 2 

=•4 



5 l 1 - *™*-!,?,^ vm)) ^ 



The result follows by (|40 |l . (|41 ]) . (|45 jl . and (j47 j) . ■ 
Theorem [Til suggests the following definition. 

Definition 7 Function Q : N — > Q n is said to oe sufficiently exciting if there 
exists a sequence of pairs of positive real numbers (e^, Ni)^L x satisfying 

/ki+Ni-l \ 

Ointa Qk)>^ ( 4 8) 

V k=ki ) 

for ki — X^7=i ^Vj with k\ — 0, and 



N 

lira Vi(£„ JVi) =oo (49) 

TV— >oo ^ — ' 



i=l 

where S(-, •) is as defined in (|40[) . 

The following result is the discrete-time analogue of Theorem @] The proof 
would have been similar to that of Theorem 31 had it not been absent from the 
paper. 

Theorem 12 Let Q : N — > Q n be sufficiently exciting. Then, for all connected 
interconnection A 6 K pxp ; solutions £j(-) of array ([55)) synchronize to 

x{k) = (r T ®/„)x(0) 

where r € K p satisfies ([4|. 

Notion of persistence of excitation carries readily to discrete time. See the 
below definition. 

Definition 8 Map Q : N — ► <2„ is said to &e persistently exciting if there exists 
a pair (e, N), e > and iV £ N>i, suc/i £/ia£ 

/ko+N-l \ 

\ fe=feo / 

/or aiZ fc S N. 
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The following theorem is the discrete-time analogue of Theorem [5] We omit 
the proof. 

Theorem 13 Let interconnection A 6 M. pxp be connected and function Q : 
N — > Q ra persistently exciting. Then solutions Xj(-) o/ array ([36]) exponentially 
synchronize to 

x(k) = (r T ®/„)x(0) 

where rel p satisfies ([4]). 

6.2 Negative results in discrete time 

Negative results generated in Subsection 14.31 are not peculiar to continuous- 
time arrays. Counterexamples similar to the ones constructed in the proofs of 
Theorem |5] and Theorem [7] can be obtained in discrete time. We thus have the 
following two theorems. 

Theorem 14 There exist maps A : N — » R pxp , where A^ is an interconnection 
for each k, and Q : N — > Q n smc/i £/iai system x + = (7 np + (A& — J p ) ® Qfc)x /ias 
an unbounded solution. 

Theorem 15 There exist connected interconnection A S M pxp and and map 
Q : N — > Q n satisfying 



such that solutions aJj(-) o/ array (|36[) do not synchronize. 

6.3 Observability grammian and synchronizability in dis- 
crete time 

As is the case with continuous-time arrays, there is a close relation between the 
observability grammian and synchronizability in discrete time. In this subsec- 
tion we will provide definitions and theorems through which we formalize that 
relation. 

Definition 9 For A : N -> R nx ™ and C : N -> R mxn , pair (C, A) is said 
to be asymptotically observable if the summand of the observability grammian, 
k i — ► ^^(k, 0)C T (k)C(k)<f>A(k, 0), is sufficiently exciting. 

Definition 10 For A : N -> R" xn and C : N -> R mxr \ pair (C, A) is said to 
6e uniformly observable i/ t/iere exists a pair (e, iV), e > and A~ G N>i, suc/i 
i/ia£ 




(W (k, k + N))>e 



for all k e N. 



2N 



The below result follows from Theorem [T2"l 

Theorem 16 Let A : N -+ R nxn and C : N -> K mxn safest Assumption^ 
(with constants a and c.) If pair (C, A) is asymptotically observable, then it is 
synchronizable. In particular, if we choose L : N — > M nx ™ as 

L(fc) := (Sc) -1 ^* + 1, 0)*5(fc, 0)C T (fc) (51) 
then for each A G <7>o solutions Xj(-) o/ array (|35[) synchronize to 

x(k) := (r T ®<f> A (k, 0)) 

where rGP satisfy 

Theorem ITBI yields the following result. 

Theorem 17 Let A : N -> R" x " and C : N -> M mxn satis/y Assumption^ 
(with constants a and c.) Then pair (C, A) is synchronizable if it is uniformly 
observable. In particular, if we choose L : N — ► tt™ xn as m (|51|) t/ien for each 
A € <7>o solutions Xi(-) of array ([35]) exponentially synchronize to 

x(k) := (r T ®<S> A (k, 0)) 

where r£P satisfy 

7 Conclusion 

We studied synchronization of stable, linear time- varying systems that are cou- 
pled via their outputs. We provided sufficient conditions on observability gram- 
mian for the existence of a bounded linear feedback law under which the systems 
synchronize for all fixed connected interconnections. Related to the main prob- 
lem, we also studied an array of coupled integrators with identical time-varying 
output matrices that are symmetric positive semi-definite. We showed, via Lya- 
punov arguments that, the trajectories of this array stay bounded. Moreover, if 
the interconnection is connected and output matrix satisfies some observability 
condition, then the systems were shown to reach consensus. 
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